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1 Introduction 

Consider a real-valued measurable zero-mean strictly stationary process Y(t), t G Z, 
obeying the following assumption. 

Assumption 1. E\Y{t)\^ < oo for all k, and Y{t) has (cumulant) spectral densities 
of orders A: = 2, 3,..., that is, there exist the functions /fc(Ai,..., Afc_i) G 
A = (—TT, tt], fc = 2, 3,..., such that the cumulant function of order k is given by 


L 


/fc(Ai,..., Afc_i)e*^i dXi... dXk-i- 


J'fc (fl T ■ • ■ : Ik—l ) 


Ak-1 


Suppose that we are given the observations where ifr = {—T ,..., 

T},T G Z. 

In this paper, we will study large sample properties of the empirical spectral func¬ 
tionals of the form 



( 1 ) 


for appropriate functions (p{X) with £ Li{A), where /^(A) is the fcth 

power of the periodogram based on the tapered data {hT{t)Y(t),t G Kt}, and k is 
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a positive integer. The taper will be of the form hT{t) = h(y) with h satisfying the 
following standard assumption. 

Assumption H. The function h{t),t € R, is an even positive function of bounded 
variation with bounded support: h{t) = 0 for |f| > 1. 

The periodogram corresponding to the tapered data is dehned as 


/t(A) 


1 


dr (A) 


Xg A, 


where dr (A) is the finite Fourier transform based on tapered data: 


driX) = 4(A) = Y. A e A, 

tGKT 


Hk,TW= [ 

J Kt 

and we suppose that id 2 ,T( 0 ) ^ 0. 

Functionals of the form (1) for fc = 1 have been extensively studied in the liter¬ 
ature, in particular, due to their applications for parameter estimation in the spectral 
domain: their behavior as T —>• oo is important for establishing asymptotic properties 
of so-called minimum contrast estimators such as Whittle and Ibragimov estimators 
(see, e.g., [9, 1], and references therein). The case of the squared periodogram was 
treated, for example, in [8], with application to a goodness-of-fit statistics, and in 
[12], with application to minimum contrast estimation. 

Asymptotic results for the functionals of the form (1) with general k > 2 were 
studied in [6] and applied to derivation of properties of weighted least squares esti¬ 
mators in the frequency domain and also in [11], with several applications discussed, 
in particular, a frequency domain goodness-of-fit testing. 

In this paper, we derive asymptotic results for functionals (1) with general k >2 
in a more general setting, using the tapered data, and under a different set of condi¬ 
tions in comparison with those used in [6] and [11]; in the Gaussian case, we state 
our results in terms of integrability conditions for the spectral density and weight 
function. Methods for the proofs are similar to those used in [12] with appropriate 
modihcations required for the more general case under consideration in the present 
paper. 


2 Results and discussion 

We begin with the following assumptions. 

Assumption!. The spectral densities fk{Xi,..., Xk-i), fc = 2, 3,..., of the stochas¬ 
tic process Y(t) are bounded and continuous. 

Assumption 3. The weight function i^(A) is bounded and continuous. 

In what follows, we denote the second-order spectral density / 2 (A) simply by 
/(A) omitting the subscript 2. 
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Theorem 1. Let Assumptions 1, 2, and H hold, and let the functions tp, (^i(A),.... 
(pk (A) satisfy Assumption 3. Then, as T ^ oo, 

( 1 ) 


EJk,T{p)^kl [ p{X)f{X)dX; 

JA 


( 2 ) COv[Jk,T{pl),Jl,T{f^2)) 

27r 


T 


e{h)klk\l\ 


J A 


(^i(A)[^(A)+^(-A)]/'=+'(A)dA 


where 


(3) 


[ V5i(Ai)</J2(A2)/^ ^(A2)/4(Ai, —Ai, A 2 ) dAi dA2 

e(^) = |y h^it)dt^ J h'^{t)dt; 

CUm{jrm,T{Pl), • ■ ■ , Jmk,T{f>k)) = 0{T^~’^). 


Suppose now that the process Y (t) is Gaussian. In this case the above asymptotic 
results can be stated under the conditions of integrability on the weight function and 
spectral density. 

Theorem 2. Let Y(t), t € Z, be a Gaussian stationary process with spectral density 
fix), X G A, such that /(A) S LpiA), and let the functions ... ,pk G LqiA), 
where 1 < p,q < 00 . Suppose also that Assumption H holds. 

(1) Ifp and q satisfy the relation 


- + k- = l, 

q p 


then, as T ^ 00 , 


EJk,Tip)^k\ [ piX)fiX)dX. 

JA 


(2) Ifp and q satisfy the relation 


1 k + l 1 _ 1 


q 


then, as T ^ 00 , 


cov[Jk,ripi), Ji,TiP 2 )) 

2'jr f 

-eih)klk\l\ / (/7i(A)[^(A)+^(-A)]/''+*(A) dA. 


T 


(3) Ifp and q satisfy the relation 


1,1 1 
—^ k- — 

q p 2 


then the cumulants of orders r > 3 of the nonnalized functionals Jk,Tif’i) tsnd to 
zero as T ^ 00 : 

cum(T^/^Jk,riPi), ■■■, T^^^Jk.riPr)) 0 . 
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(4) Ifp and q satisfy the relation 

1 ^ ki ‘ k>p 1 1 

^ {*) 5 

q r P ^ 

then the cumulant of rth order, r > 3, of the normalized functionals Jki,T{pi), 'I = 
1,..., r, tends to zero as T ^ oo: 

cum{T^^'^Jki,T{<Pi), ■ ■ . ,T^^^Jkr,Ti‘Pr)) 0. 

As corollaries of the above theorems, we obtain the next asymptotic normality 
results. 

Let us fix the weight functions ipi,..., ipm and denote 
Jt = and J = 

where Jk{ip) = k\ ip(X)f^(X) dX. 

Let ^ be a Gaussian random vector with zero mean and second- 

order moments 

Wij = E^i^j 

= 2TTe(h)(kky (pi{X)[ifj{X) + ifj{-X)]p'‘{X) dX 

+ J <Pi(Ai)<Pj(A2)/^ ^(Ai)/^ ^(A2)/4 (Ai, —Ai, A 2 ) dAi dA2^ . 

Assumption 4. The spectral density of the second order /(A), the weight function 
(/3(A), and the taper h are such that T^/'^{EJk^T{v) ~ Jk{^)) 0. 

Theorem 3. Let Assumptions 1, 2, and H hold, and let the functions ipi, i = 1,... ,m, 
satisfy Assumption 3. Then 

- EJt) ^ ^ asT^oo; 

moreover, if Assumption 4 holds for the functions ipi, i = 1,... ,m, then 

asT^oo. 

Let (( = {Ci}i=i....,m be a Gaussian random vector with zero mean and second- 
order moments 

Vij = EQ^j = 2TTe{h){kk\Y f (p,{X)[(pj{X) +‘Pj{-X)]y'^{X) dX. 

J A 

Theorem 4. Let Y(f), t (2 Z, be a Gaussian stationary process with spectral density 
/(A) G Lp{A), and let the functions pi ,..., pm G Lq(A), where 1 < p,q < oo, be 
such that - kX = X Suppose also that Assumption H holds. Then 

q p z 

- EJt) ^ C asT^oo; 

moreover, if Assumption 4 holds for the functions pi, i = 1,... ,m, then 

J) ^ C asT^oo. 
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We next state more general results, the joint asymptotic normality for functionals 
of different powers of the periodogram. 

Consider 

Jki,...,k^,T = and Jki,...,krn = 

where Jki^p) = k\ ip{X)f'^{X) dX. 

Let ^ be a Gaussian random vector with zero mean and second- 

order moments 

Wij = 

= 2Tre{h)kih\kjkjl^J [(^j(A) + (pj{-X)]p^^'"^ {X) dX 

+ [ ^^iXl)^J{X2)f^-\Xl)P^-\X2)f4iXu-XlA2)dXldX2 
JA'^ 

Theorem 5. Let Assumptions 1, 2, and H hold, and let the functions ipi, i = 1,... ,m, 
satisfy Assumption 3. Then 

T^^^{Jki,...,km,T — EJki,....km,T) ? asT ^ oo; 
moreover, if Assumption 4 holds with k = ki, ip = ipi, i = . ,m, then 

T^^‘^iJki,...,k,„,T - Jki,...,k,„) ^ C asT ^ oo. 

Let = {Ci}i=i....,m be a Gaussian random vector with zero mean and second- 
order moments 

= ECf = 2TTe{h)hh\kjkf. f pi{X) [^.(A) + ^j(-A)] (A) dX. 

J A 

Theorem 6. Let Y(t), t G Z, be a Gaussian stationary process with a spectral 
density /(A) G Lp{A), and let the functions pi,..., pm G Lq{A), where 1 < p, g < 
oo, be such that A -f minj/cijA = A. Suppose also that Assumption H holds. Then 

T^^'^{Jki,...,k,„,T — EJki....,k,n,T) C asT ^ oo; 
moreover, if Assumption 4 holds with k = ki, p = pi, i = 1,... ,m, then 

T^^‘^{Jki,...,km,T — Jki,...,km) C as T ^ oo. 

Remark 1. Integrals of nonlinear functions of the periodogram (including, in par¬ 
ticular, powers of positive orders of the periodogram) were studied, for example, in 
[13] for discrete time processes under the assumption of boundedness of the spec¬ 
tral density. In [8], the integral functionals of the squared periodogram were studied 
for stationary Gaussian series given by the moving-average representation, and the 
asymptotic normality result was stated under the particular assumption of summa- 
bility of the coefficients of the representation and continuity of the derivative of the 
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spectral density. In [6] and [11], the asymptotic results for functionals of powers of 
the periodogram of general order have been studied under the conditions of summa- 
bility of cumulants of the process. In this paper, we state the results for discrete-time 
non-Gaussian processes under the condition of boundedness of spectral densities of 
all orders (which are supposed to exist), and we also derive the results for Gaussian 
case under the conditions of integrability of the spectral density and weight func¬ 
tion. 

Remark 2. Conditions on the spectral density under which Assumption 4 will be 
satisfied can be formulated analogously to the corresponding conditions in [1] for the 
case where h{t) = 1 and analogously to the conditions in [2] for the general h{t) of 
Assumption H. 

Remark 3. The results on asymptotic properties of the integrals of the powers of the 
periodogram can be useful for some problems of statistical inference. One possible 
application is hypothesis testing concerning the form of the spectral density of the 
process. For example, in [8], a quadratic goodness-of-fit test in the spectral domain 
was studied for Gaussian processes. Note that the asymptotic normality result for the 
corresponding test statistic stated in [8] can be also derived from our Theorem 6, 
that is, under a different set of conditions. More applications of the integrals of the 
powers of the periodogram for goodness-of-fit testing, peak testing, and assessing 
model misspecification are presented in [11]. In [6] and [12], the integrals of the 
squared periodogram were applied for parametric estimation in the spectral domain. 

3 Proofs 

For the proofs of the results of Section 2, we use the technique based the properties 
of the multidimensional kernels of Fejer type (see, e.g., [5, 1], and references therein) 
and the Holder-Young-Brascamp-Lieb inequality (see [3, 4], and references therein; 
see also [12]). In what follows, we will refer to the latter as the HYBL inequality. The 
application of these tools leads to very transparent and elegant proofs. We will also 
use the formula giving expressions for cumulants of products of random variables via 
products of cumulants of the individual variables (see, e.g., [10, 5]) and the multi¬ 
linearity property of cumulants. The lines of reasonings are very close to those used 
in [12]. 

Proof of Theorem 1. Consider 

^ [icum{dTiX)dT{-X)) ^ 

= ■ 

We apply now the formula for cumulants of products of random variables (see, e.g., 
[10]); it is convenient to assign the indices to As in the following way: we can enu¬ 
merate all As appearing in the above row from 1 to 2k, having in mind that Ai with 
odd indices are simply equal to A, whereas Ai with even i are equal to — A. Then we 
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can write down the expectation in the following form: 

p 


e4{x) = 


1 


(27r7?2,T(0))'= 


E n cum{dT{\i),i G vi) 


partition of (1,...,2A:) 


S{X2i-i — A) <5(A2i + A). 


( 2 ) 


The cumulants of the hnite Fourier transforms dr (A), X G A, can be written as 
follows: 

cum{dT{ai),.. .,dT{ak)) 

P k 

= / (.k{ti-tk,...,tk-i-tk)dti...dtk 

r=i 

= / /fc(7i,---,7fc-i) 


lAk-l 

k-1 


k-1 


ddi,r(7i “E^'i ““fe d'ji.. .djk-i, 


i=i 


i=i 


where 

Hi^t{X)= [ hT{t)e-^*^ dt. 

J K'j' 

Correspondingly, we obtain the following formula for the expectation of Jk^Ti^p)- 


EJk,T{p)=E [ p{X)4iX)dX 
Ja 

= [ pW 

Ja 


(27r7T2,T(0))' 


E 


partition of (l,...,2fc) 


f 


YlfWiliijJ e i>z) 


2k 


P 


X n ^i.t(7j - n ( e ) n <5(A2^-l - A) n + A) dy dX. 

j=l 1=1 ' i=l i=l 

(3) 


Here and in similar formulas below, we use the following notation: for a set of 
natural numbers v, we denote by \v\ the number of elements in v and by v the subset 
of V that contains all elements of v except the last one. Integration in the inner integral 
in the above formula is understood with respect to {2k — p)-dimensional vector 7 ' 
obtained from the vector 7 = ( 71 ,..., 72 ^) due to p restrictions on the variables 7 ^, 
j = 1,..., 2k, described by the Kronecker delta functions 5. 

Now we note that the products 0^=1 ^i,t(Aj) in the case where 'Yf^i = 0 
give rise to a class of d-type kernels (or Fejer-type kernels). Namely, if Assumption H 
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holds and i?fe,T(0) ^ 0, then 




(2^)fe-ii7fe,T(0) 


k-l 

i=i 


' fc -1 

^ i=i 


(4) 


is a kernel over which is an approximate identity for convolution (see, e.g., 

[7]), and 


lim / G'(ui 
T^ooJAk-^ 

= G{vi ,... ,w/c-i), 


(5) 


provided that the function Gis bounded and continuous at the point {vi,, 
Vk-l)- 

The asymptotic behavior of the right-hand side of (3) can be evaluated basing on 
the property (5). 

Let us first consider the partitions ly composed by pairs. For those partitions, when 
the products of only the cumulants of the form cum{dT{X), dT{—X)) appear in (2), 
we obtain under the integral sign in (3) the terms of the form 

= {y /(7)^E(7- A)rf7| ■ 


We note that there are fc! such terms, therefore, in the expression for EJk^xi^), we 
have the term ^ 

fciy v?(A)|y/(7)^>^ ^(7-A)d7} d\ (6) 

and this is the only case where we have k kernels, and all k factors 2-nHl^T(o) 
to compose these kernels ^2 t(')- 

In all other partitions, we will be able to compose from 1 to fc — 1 kernels taking 
combination of with suitable arguments: for each 2nd-order kernel, we will 

use one of the factors 2 -khIt(o) ^om ^27177T(o)T^ ’ otherwise, when composing the 
hh order kernel with Z ^ 2, we will need the normalizing factor ^ 2 -Ky-^Hi t(o) ’ 

therefore we will modify the factor ( 2 ^g,^^(o))'= ^y taking, instead, ■ 

So, for those partitions, when we compose kernels of orders, say, /i,..., 4, with 
Si=i h = 2fc, the corresponding integral in (3) will be represented in the form of a 
generalized convolution of some product of spectral densities of different orders with 
the product of kernels of orders li,... Jr, and the factor 

(27riL2,T(0))'= 


( 7 ) 
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will be supplied to the integral. For example, in the simplest case where p = \, the 
corresponding term in (3) can be represented as follows; 


( 2 ^)^fc-iH 2 fc.T( 0 ) 

( 2 ^i/ 2 ,T( 0 ))'= 

X ^2fc,T(7l ~ ^ 

k 

X l^(A 2 i-l — 


ip{X) 


! J\2k-1 

1 ) ■ • ■ , l2k-l — X2k-l) 


/2fe(7l,---,7fc-l) 


k 2k-l 

A) n +A) n dx. 

i—l i—1 


Now we take into account the following asymptotics for Hk^riO)' Hk^riO) ^ 
THk{0), where Hk{0) = f /i^(A) dX, and conclude that, in the case of r kernels, 
1 < r < fc — 1, the factor (7) is asymptotically of order ; the corresponding 
integrals containing these kernels will converge to hnite limits under the conditions 
of the theorem according to (5). Therefore, the expectation is obtained as the limit 
of ( 6 ). This gives statement (1) of the theorem. 

Consider 


COv[jk,T{^Pl), Ji,t{V2)) 

1 

“ (27ri72.T(0))'=+' 

X [ (pi{a)T^i(3)cum{{dTia)dTi-a)f,{dTiP)dTi-P))'')dadp. (8) 

Ja^ 

The cumulant under the integral sign in ( 8 ) according to the formula for calculation 
of cumulants of products of random variables can be written in the form 

p 

E n cum{dT{p.j),Pj &Vi), (9) 

v={v-i_,...,Vp) i=l 

where the summation is taken over all indecomposable partitions v = (z/i,..., Vp), 
|z/i| > 1 , of the table T 2 with two rows, {a, —a ,... ,a, —a} (of length 2k) and 
{/3, —/3,..., /3, —/3} (of length 21). For asymptotic analysis of expression ( 8 ), we can 
use the reasonings analogous to those for the case of functionals of squared peri¬ 
odogram in [ 12 ], but now, dealing with the tapered case, we need to keep track of 
normalizing factors for appearing kernels. Again, similarly to the previous consid¬ 
eration of the expectation, we analyze all possible partitions and kernels that can be 
composed for every particular partition. Let us hrst consider the terms in (9) that 
correspond to partitions by pairs, that is, 

k+l 

Y\_cum{dT{pi),dTiX^)), (10) 

i=l 

where pi, Xi G {a, —a, j3, —/3}, and v = {(pi, Ai), i = 1,... ,k + 1} forms an inde¬ 
composable partition of the table T 2 . 

In the case where we have the k — 1 cumulants cum{dT{a), dri—a)) and I — 1 
cumulants cum{dT{l3),dT{—0)) in the product (10), according to formula (4), we 
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can compose k + I — 1 kernels {k + I — 2 kernels of order 2 and one of order 4), and 
the factor before the integral in ( 8 ) becomes of the form which is asymp¬ 

totically of order y. Only the terms of this kind in (9) give the main contribution (of 
order into the covariance ( 8 ); all other terms in ( 10 ) produce a smaller-order con¬ 
tribution to the covariance ( 8 ). More precisely, in order to describe the asymptotics of 
the covariance, we have to consider, among the terms in (9), the following ones; 

(cum(dT(a), dr (—a))') ^ ^ (cum(dT(l3),dT(-/3))y ^ 

X ^cum(dT(o:), dT(/ 3 )) cum(^dT (-a),dT(-/ 3 )} 

+ cum (dr (a), dr (—/3)) cum (dr (-a),dT(/3)}]. (11) 


Their contribution to the covariance is of the form 

J /(7i)iTi,T(7i - a)-f^i.T(-7i + a) c^7i 

J f { 72 )Hi{'12 - P)Hi^t{- 12 + P) di2 
J /(73)-ffi.T(73 - a)Hi^T{-73 - P) ^73 

X j /(74)-ffi,T(74 + a)-ffi,r(-74 + ;3) ^74 

+ J f{l3)Hl,T{73 -ol)Hi^t{- 73 +P)di3 

X j f{74)Hi^T{74 + a)Hi^T{-74-P)di4 


-I k-l 


1 i-1 


da dp 


27riJ4,T(0) 

{H 2 ,t{ 0 W 


J J ‘Pi{a)‘P2{P) y /(7i)^2,t(7i-Q;) d7i 


k-l 


n 1-1 


/(72)<^2,t(72 - ^) ^72 


'4.t(73 - a, -73 - P, 74 + a) 


+ ^4 ,t(73 - a , -73 + P,"f 4 + a)] ^73 ^74 da dp. 

Taking into account formula (5), we can evaluate the latter as 

We note also that there are klkW. terms of the form (11) in (9). 

Let us now consider the terms in (9) that correspond to other partitions. We can 
see that one more possibility is left to compose k + l — 1 kernels, namely, the case of 
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the terms of the form 

[cum(dr(a), dra)) f ^ [cum{dTiP),dTi-/3))]’' ^ 

X cum{dT { a ), dT {- a ), dTi / 3 ), dri - jd )) 

in the sum (9) (there are klkUl such terms), and the corresponding contribution to the 
covariance (8) is of the form 


X y/(7i)-ffi.T(7i - a)-f^i.T(-7i + a) c^7i 


-I k-l 


1 i -1 


J /(72)-ffl,T(72 - ,S)iJl,T(-72 + ,3) ^72 

X Hi^rips - /3)^f X! ^r-2 dfj ,3 da dfd 

27riJ4.T(0) f f , 

y ^4(aV2(/3) 

X J /(7i)^'2,t(7i- a)(i7i J /(72)^2.r(72 -/3) ^72 
^ 11J Ct,d 2 +a, fj,3-13) d^id^i 2 dfi 3 dadfd 

Y (J fe2(yyy J J Vi{a)W2{P)f{a)f{P)fi{oi,-a,P)dadP asT^cx). 


In all other cases, we can compose less than k + l — 1 kernels; the corresponding inte¬ 
grals will converge to hnite limits supplied by the factor of orders not exceeding ^. 

Summarizing the above reasonings, we come to the asymptotics for the covari¬ 
ance as given in statement (2) of the theorem. 

We now evaluate the asymptotic behavior of the cumulant of order k > 3: 


cum{j^^^Tiv>i) , ■ ■ ■ , Jrrik ,T iV>k)) 

^ (2^172,t( 0 ))^ 

X cM77i((dT(ai)dT(-ai))™\ . • ■, {dT{ak)dT{-ak))""°) dai . ..dak, 
where M = Yf=i 

The cumulant under the integral sign can be represented as the sum 


p 

E n cum[dT{lJ.j), l^j G Vi), 

,^=(ui,...,Vp) i=l 


( 12 ) 
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where the summation now is taken over all indecomposable partitions z/ = (i^i,..., 
z/p), jz/il > 1 , of the table Tk with k rows {ai, —at, ..., Ui, —ai}, i = 1 ,... ,k, 
the length of the zth row being 2mi. Starting again with consideration of partitions by 
pairs, we can see that with these partitions we can compose at most ~ 

1) + 1 = M — k + 1 kernels (M — k kernels ^ 2 T kernel 

the corresponding integrals will converge to finite limits supplied with the factor 
= ( 2 rS’( 0 )°^ ’ '^hich is asymptotically of order With 

all other partitions, we will be able to compose no more that M — fc + 1 kernels; 
therefore, their contribution to cumulant (12) will be of order less than This 

gives statement (3) of the theorem. □ 

Proof of Theorem 2. We use the same calculations as those in the proof of Theo¬ 
rem 1 , but to analyze the limit behavior of the integrals representing the cumulants, 
we will appeal to the HYBL inequality (see [3, 4, 12]). The reasonings follow the 
same lines as in [12], so here we just point out the key steps. 

For the Gaussian case, we will have only partitions by pairs in (3): 


EJkM = E f ip{X)I^iX)dX 

JA 

= f v’W 

JA 


^(27riT2,T(0))'= 




partition of ,2k) 


2k 


k 


X n n ) n - a) n <^(A2 *+a) dy dx. 

j=l Z=1 Ygi/, / j—1 i—i 

(13) 


We consider separately the term ( 6 ): 


fc!^V5(A)|^/(7)<^2.T(7-A)d7} dX 


= k\ 


! A^ 


‘F(A) Jl/(Ti - A) dX 


djj. 

i=i i=i 


Note that the convergence to the finite limit k\ ip(X)f^(X) dX will be assured if we 
assume the conditions for statement ( 1 ) of the theorem. 

Now consider the remaining terms: we have the integrals over with inte¬ 
grands composed by products of the functions (p with k functions / and 2 k functions 
Hi^t with some linear relations between the arguments of these functions; these in¬ 
tegrals are supplied with the factor ( 27 rg 2 !r(o))’^ ' 

Applying the HYBL inequality, we can bound each such integral by the expres¬ 
sion 


1 


{ 2 nH 2 .Tm 




( 14 ) 











Asymptotics for functionals of powers of a periodogram 


193 


provided that (p{X) G Lq{A), /(A) G Lp(A), and iJi t(A) G Lr(A) with 

— -j- k — -t- ‘2k — = /c “t" 1. 

q p r 

If we choose r = 2 and take into account that, under Assumption H, we have 
||-ffi,T||r < and ~ T, then from (14) we arrive at the bound 

consfllt/jllqll/llp as T ^ c», with the restrictions onp and q as in statement (1) of the 
theorem. From this point we can repeat the same arguments as in [12] to show that, 
in fact, this bound can be strengthen to o(l) as T —> oo. 

The similar reasonings are applied to derive statements (2)-(4) of Theorem 2. □ 
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